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Abstract 

It is well known that in fusion reactions one may get only neutron deficient superheavy nuclei 
located far from the island of stability. The multi-nucleon transfer reactions allow one to produce 
more neutron enriched new heavy nuclei but the corresponding cross sections are rather low. 
Neutron capture process is considered here as alternative method for production of long-lived 
neutron rich superheavy nuclei. Strong neutron fluxes might be provided by nuclear reactors and 
nuclear explosions in laboratory frame and by supernova explosions in nature. All these cases are 
discussed in the paper. 

PACS numbers: 25.70.Jj, 25.70.Lni 
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I. MOTIVATION 



A ten years epoch of ^^Ca irradiation of actinide targets for the synthesis of superheavy 
(SH) elements is over. The heaviest available targets of Berkelium [Z = 97) and Californium 
{Z = 98) have been used to produce the elements 117 Q and 118 Q- Note that the predicted 
cross sections and excitation functions j^, 4| have been fully confirmed by the experiments 
performed in Dubna and later in Berkeley and GSI. These predictions are based on the shell 
corrections to the ground states of superheavy nuclei calculated by P. MoUer and others 
jsl which we used to estimate survival probabilities of these nuclei against fission. It is 
the gradual increase of the fission barriers of the compound nuclei (CN) formed in these 
reactions which explains the almost constant value (of a few picobarns) of the cross sections 
for the production of SH elements with Z = 112 118 in hot fusion reactions. 

As mentioned above, 249Qf (T1/2 = 351 yr) is the heaviest available target that can 
be used in experiment. Thus, to get SH elements with Z > 118 in fusion reactions, one 
should proceed to heavier than ^^Ca projectiles. The strong dependence of the calculated 
evaporation residue (EvR) cross sections for the production of element 120 on the mass 
asymmetry in the entrance channel {g] makes the nearest to ^^Ca projectile, ^'^Ti, most 
promising for further synthesis of SH nuclei. Of course, the use of the titanium beam instead 
of ^^Ca also decreases the yield of SH nuclei mainly due to a worse fusion probability. The 



calculated excitation functions for the synthesis of SH elements 119 and 120 in t 
reactions of ^°Ti with ^^^Bk and ^^^Cf targets are shown in Fig. [1] taken from Ref. 
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FIG. 1: Excitation functions for production of superheavy elements 119 and 120 in the 3n and 4n 
evaporation channels of the ^'^Ti+^^^Bk, ^"^Ti+^^^Cf and ^^Cr+^^^Cm fusion reactions. 



The orientation effects are known to play an important role in fusion reactions of statically 
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deformed heavy nuclei [sl,^]. The fusion probabihty (formation of CN) is strongly suppressed 
for more elongated nose-to- nose initial orientations. As a result, the preferable beam energies 
for the synthesis of SH elements in the hot fusion reactions are shifted to values that are 
several MeV higher than the corresponding Bass barriers (calculated for spherical nuclei). 
As can be seen from Fig. [H the estimated EvR cross sections for the 119 and 120 SH elements 
synthesized in the ^°Ti induced fusion reactions are quite reachable at available experimental 
setups, though one needs much longer time of irradiation than for the ^^Ca fusion reactions. 

The yield of superheavy nuclei (number of events per day) depends not only on the cross 
section but also on the beam intensity and target thickness. In this connection the other 
projectile-target combinations should be also considered. Most neutron-rich isotopes of ele- 
ment 120 may be synthesized in the three different fusion reactions ^^Cr-|-^^^Cm, ^'^Fe+^^^Pu 
and ^^Ni+^'^^U leading to the same SH nucleus ^°^120 with neutron number near to the pre- 
dicted closed shell N = 184 [3]. These three combinations are not of equal value. The 
estimated EvR cross sections for more symmetric ^^Fe+^'^^Pu and ^^Ni+^^^U reactions are 
lower than those of the less symmetric ^^Cr+^^^Cm combination js], which, in its turn, is 
quite comparable with the Ti-induced fusion reaction, see the right panel of Fig. [H The 
advantage factor 2 or 3 for the ^°Ti+^^^Bk and ^°Ti+^^^Cf fusion reaction as compared with 
^^Cr+^^^Cm is definitely within the theoretical uncertainty for calculation of such small cross 
sections. 

Such uncertainty originates from two main factors. First, it is connected with the complex 
dynamics of CN formation in competition with the dominating quasi-fission process. There 
are too few experimental data on low-energy collisions of very heavy ions to fix finally all 
the parameters (potential energy, friction forces, nucleon transfer rate) and to obtain an 
explicit value for the probability of CN formation for such heavy systems. The second 
source of uncertainty comes from the estimation of the survival probability of the excited 
CN, which is calculated within the statistical model. The basic quantities here are the fission 
barriers of SH nuclei appearing in the neutron evaporation chain. For example, in the 4n 
evaporation channel one needs to know the fission barriers and neutron separation energies 
for 5 subsequent isotopes of SH element. 

Different theoretical models give quite different values for the fission barrier for heavy 
nuclei. As mentioned above, we used the fission barriers of SH nuclei predicted by the 
macro-microscopic model jsj, which gives lower fission barriers for the isotopes of elements 
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FIG. 2: Excitation functions for the production of superheavy elements 119 and 120 in the 4n 
evaporation channels of the ^'^Ti+^^^Bk, ^'^Ti+^'^^Cf and ^^Cr+^^^Cm fusion reactions. A factor 3 
is assumed for the uncertainty of the survival probability caused by the uncertainty of the fission 
barriers of SH nuclei. 

119 and 120 than for elements 114 -^ 116 synthesized in "^^Ca fusion reactions. On the 
other hand, the full microscopic models based on self-consistent Hartree-Fock calculations 
usually predict much higher fission barriers for these isotopes (up to 10 MeV) if the Skyrme 
forces are used [8|] (though these predictions are not unambiguous and depend strongly on 
chosen nucleon-nucleon forces). The exponential dependence of the survival probability on 
the values of the fission barriers leads to a change of it approximately by a factor 10 if the 
barriers are changed by 1 MeV. 

Thus, we may conclude that the theoretical uncertainty in the predictions of the cross 
sections for formation of SH nuclei in hot fusion reactions with actinide targets is not less 
than a factor 3. If we consider such an uncertainty in our calculations, the excitation 
functions for synthesis of the new SH elements 119 and 120 in 4n evaporation channels of 
the 50Ti+249Bi^^ 50rji^249Qf 54q^^248q^ f^gj^^ reactions will look as shown in Fig. H 

From this figure it is clear that all these reactions must be considered as quite promising 
for the synthesis of new elements. The final choice between them depends not so much 
on the difference in the cross sections as on other experimental conditions (availability of 
appropriate targets, beam intensities, etc.). 

Note, however, that the present limits for neutron rich isotopes located at the upper part 
of the nuclear map (Z > 60) are very close to the line of stability while the unexplored area 
of heavy neutron-rich nuclides (to the east of the stability line) is extremely important for 




FIG. 3: Upper part of the nuclear map. The island of stability is shown schematically. 

nuclear astrophysics investigations and, in particular, for the understanding of the r-process 
of astrophysical nucleogenesis. For nuclei with Z > 100 only neutron deficient isotopes 
(located to the left of the stability line) have been synthesized so far, see Fig. [31 Due to the 
bending of the stability line to the neutron axis, in fusion reactions of stable nuclei one may 
produce only proton rich isotopes of heavy elements (it's true also for the isotopes of 119 
and 120 elements discussed above). This is the main reason for the impossibility of reaching 
the center of the "island of stability" in the superheavy mass region {Z ~ 114 120 and 
~ 184) in fusion reactions with stable projectiles. In the ^'^Ti+^'^^Cf and ^^Cr+^^*Cm 



fusion reactions the 120-t 
0.02 pb, correspondingly 



1 element can be synthesized with cross sections of about 0.04 and 







but again only the short-lived neutron deficient isotopes of this 
element may be obtained. The use of beams of radioactive nuclei hardly may solve this 
problem due to their low intensities. 

Multi-nucleon transfer processes in low-energy collisions of actinide nuclei (like U+Cm) 
may really lead to the formation of neutron rich long-living superheavies. The shell effects 
(antisymmetrizing quasi-fission process with a preferable formation of nuclei close to the 
doubly magic lead isotope) might significantly enhance the corresponding cross sections 
6|, l9| . In spite of difficulties of separation of the transfer reaction products, experiments of 
such kind are planned to be performed in nearest future. 

The neutron capture process is an alternative (oldest and natural) method for the produc- 
tion of new heavy elements. Strong neutron fluxes might be provided by nuclear reactors and 
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nuclear explosions under laboratory conditions and by supernova explosions in nature. How- 
ever the "Fermium gap" , consisting of the short lived Fermium isotopes 258-260pj^ located 
at the beta stability line while have very short half-lives for spontaneous fission, impedes 
formation of nuclei with Z>100 by the weak neutron fluxes realized in existing nuclear reac- 
tors. In nuclear and supernova explosions (fast neutron capture) this gap may be bypassed, 
if the total neutron fiuence is high enough. Theoretical models predict also another region 
of short lived nuclei located at Z=106-^108 and A~270 (see Fig. H]). 




120 



115 



^1. 




150 



160 



170 



180 



FIG. 4: Half-lives (up) and decay modes (bottom) of nuclei in the upper part of the nuclear map. 
Schematic view of slow (terminated at the short-lived fission Fermium isotopes) and fast neutron 
capture processes with subsequent beta-minus decays are shown by the arrows. 

In this paper we study the possibility of synthesizing heavy elements in multiple "soft" 
nuclear explosions and in pulsed reactors. We have found that in the first case the both 
gaps may be easily bypassed and, thus, a measurable amount of the neutron rich long-lived 
superheavy nuclei of the island of stability may be synthesized. For the second case we have 
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formulated requirements for the pulsed reactors of the next generation which could be also 
used for the production of long-lived superheavy nuclei. 



II. NUCLEOSYNTHESIS BY NEUTRON CAPTURE 



The synthesis of heavier nuclei in the reaction of neutron capture with subsequent beta- 



minus decay is a well studied process (see, e.g., |10l-ll2|). Relative yields of the isotopes 
formed in such a process may be found as a solution of the following set of differential 
equations (somewhat simplified here) 



where no is the neutron flux (number of neutrons per square centimeter per second) and 
X^z A ~ ^^'^/'^i/2 decay rate of the nucleus {Z,A) into the channel i (i.e., beta-minus 

and alpha decays and fission). For simplicity we ignore here the energy distribution of 
the neutrons and, thus, the energy dependence of the neutron capture cross section af^^. 
Neutrons generated by fission in nuclear reactors and in explosions are rather fast (far from 
the resonance region). In the interval of 0.1-1 MeV the neutron capture cross section is a 
smooth function of energy with the value of about 1 barn, which is used below for numerical 
estimations. Note that integration over the neutron energy may be performed very easily 
and does not change the obtained conclusions. The simplest version of the reaction chain, 
where only the neutron capture reactions are retained [first line in Eq.(Il])], was considered 
recently in Ref . [isl . 

To solve Eq. ([1]) numerically one needs to know the decay properties of neutron rich nuclei 
which are not studied yet. That is the main problem, which significantly complicates the 
analysis of the multiple neutron capture processes. Theoretical estimations of half-lives for 
a-decay are rather reliable because they depend only on ground state masses which are very 
close in different theoretical models. For the half-lives of a-decays we used the well-known 

n 

Viola-Seaborg formula with the coefficient proposed by Sobiczewski et al. [14^ Half-lives of 
allowed /3-decays depend also on the ground state masses of the nuclei and may be estimated 
by the empirical formula log^o i^^'^) = 5-7, where the Fermi function /o is calculated 



by the standard formulas (see, e.g., Ref. |15|). 

The fission half-hves are the most uncertain quantities in Eqs.(IT]). Rehable analysis of 
the fission process requires knowledge of the multidimensional potential energy surface as 
well as of the collective inertia parameters. Such calculations can be performed only in very 
restricted areas of the nuclear map. We used the ground state masses and shell corrections 
for heavy and superheavy nuclei proposed by P. MoUer et al.[5] and then we applied the 
empirical formula for the estimation of the fission half-lives. For this purpose we employed 



the relation by Swiatecki 16| based on the idea of the dominant role of the fission barrier 
for the fission probability 

logioT/^2 (sec) = 1146.44 - 75.3153^7^ + 1.63792 {Z^/Af 
-0.0119827 {Z'^/Af + Bf (7.23613 - 0.0947022^7^) + 

0, Z and N are even 
+ < 1.53897, A is odd (2) 
0.80822, Z and N are odd 

Here Bf = B^^'^^ + 5Ug,s. is the fission barrier, which was calculated as a sum of the liquid- 
drop barrier and the ground-state shell correction. The coefficients of the systematics ([2]) 
were determine d by the fitting to the experimental data and to the rather realistic theoretical 
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18| for the region of 100 < Z < 120 and 140 < < 190. 



predictions 

The results of such calculations are shown in the left part of Fig. HI where the theoretical 
values of ^^"2' "^1/2 ^i-nd/or Tj^^ were replaced by experimental ones if known. In accordance 
with our calculations, the most stable superheavy nuclei located at the island of stability 
(which could be searched for in nature) are the isotopes of Darmstadtium {Z = 110) and/or 
Copernicium {Z = 112) with the neutron number ~ 180. To produce these nuclei in 
the neutron capture process one needs to bypass the two area of short-living fissile nuclei, 
namely, the Fermium gap {Z = 100) and the region of Z=106-^108 and A~270. 

To test our model we first described available data on the fast neutron capture process 
realized in nuclear explosions. In this case the time of neutron capture, r„ = (noa^-^)"^ ~ 
1 /is << Ti/2{Z, A), is much shorter than half-lives of the produced nuclei (up to the neutron 
drip line). Keeping only the first two terms of the r.h.s. of Eq.([T]), we get the following 
analytical solution (with initial conditions Nz,A{t = 0) = 1 and Nz,A+k(t = 0) = at A; > 
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where k is the number of captured neutrons) 

Nz,A+k = 



X 

IT 



(3) 



This relation can be used for clear understanding of the fast neutron capture process and 
for a preliminary estimation of relative yields of heavy nuclei synthesized in such a process. 



Here x = na. 



n-yi 



n 



riQT is the total neutron fluence (neutrons per square centimeter) and r 



is the duration of explosive neutron irradiation. Thus, the dimensionless quantity x = na-n^y 
is a key factor characterizing the neutron capture process. In nuclear explosions the neutron 
fluence reaches the values of 10^^ cm~^, so that a; ~ 1 and more than 10 neutrons could be 



captured during one exposure with time duration of about 1 /xs 




FIG. 5: Experimental (open dotes) and calculated relative yields of heavy nuclei in the test nuclear 
explosion "Mike"Q. 

In Fig. [5] the experimental data on the yield of transuranium nuclei in the test thermonu- 
clear explosion "Mike" {l9j] are compared with those calculated by Eqs.(IT|) assuming 1 fis 
neutron exposure of 1.3 x 10^'^ neutrons/cm^ with subsequent one-month decay time. Note 
that elements 99 and 100 (Einsteinium and Fermium) were first discovered just in debris of 
the "Mike" explosion. As can be seen, in this case the Fermium gap does not influence the 
yields of nuclei with Z > 100 which follow roughly the relation (j3]). 

Recently we proposed to consider the possibility of generating two or several nuclear 
explosions in close proximity of each other to increase the resulting mass number of the 



synthesized nuclei ISj. Here we study for the first time such a possibility illustrated in the 
upper part of Fig. [61 In the bottom part of this figure the probabilities of heavy element 
formation are shown for one, three and ten subsequent short-time (1 /xs) neutron exposures 



of 10^^ n/cm^ each following one after another with time interval of 10 seconds with final 
one month waiting (needed to perform some experimental measurements). 

Our results demonstrate for the first time that multiple rather "soft" nuclear explosions 
could be really used for the production of noticeable (macroscopic) amount of neutron rich 
long-lived superheavy nuclei. Leaving aside any discussions on possibilities of such processes 
and associated technical problems, we want to emphasize a sharp increase of the probability 
for formation of heavy elements with Z > 110 in the multiple neutron irradiations (enhance- 
ment by several tens of orders of magnitude, see Fig. |6]). This probability is high enough for 
some superheavy elements (see the region above the dotted line in Fig. E]) to perform their 
experimental identification. 




160 180 
number of neutrons 





10x(l |is + 10s) 
+ 1 month decay 


\ 1 
\ 1 
t 1 
t 

«-« 

single explosion: 

lO^'^n/cm^ 

+ 1 month decay 


I /A 

-J I 1 
1 ( 
\ 1 

■ * ' 

\ /\ i';3x(l |is+ lOs) 
\.^t 1 • *i + 1 month 

. . ...... 



100 105 110 

charge number 



FIG. 6: Schematic picture for multiple neutron irradiation of initial ^^^U material (up) and prob- 
ability for formation of heavy nuclei (bottom) in such process (one, three and ten subsequent 
explosions). Dotted line denotes the level of few atoms. 

It would be interesting also to study the same process of multiple neutron exposures 
realized in pulsed nuclear reactors. The pulse duration here could be much longer than in 
nuclear explosions (up to few milliseconds). In spite of that, the neutron fiuence usually 
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does not exceed 10^^ n/cm^ in existing nuclear reactors (no ~ 10^^ n/cm^s during one 
millisecond pulse). Thus, the quantity x = nan-y is about 10~^, the time of neutron capture 
Tn = (noan-y)^^ ~ 10^ s, and only the nearest long lived isotopes (A+1 or A+2) of irradiated 
elements may be formed during the pulse, see formula Multi-pulse irradiation here 
corresponds, in fact, to the "slow" neutron capture process, in which new elements with 
larger charge numbers are situated close to the line of stability and finally reach the Fermium 
gap where the process stops (see Fig. H]). The result of numerical solution of Eq.(IT]) for 
neutron capture process in an ordinary pulsed reactor is shown in Fig. [7] by the dashed 
line. The probability for formation of heavy elements with Z > 100 is negligibly small here 
independent of a number of pulses and total time of irradiation. 
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FIG. 7: Relative yields of heavy elements in ordinary (dashed line) and in high-intensity pulsed 
reactors (solid line) at the same total neutron fluence 10^^ n/cm^. 

The situation may change if one would be able to increase somehow the intensity of the 
pulsed reactor. The neutron fluence in one pulse and frequency of pulses should be high 
enough to bypass the both gaps of short-lived nuclei on the way to the island of stability 
(see left part of Fig. H]). Thus, the specification of the high-intensity pulsed reactors of next 
generation depends strongly on properties of heavy neutron rich nuclei located to the right 
of these gaps. These nuclei are not discovered yet (see Fig. |3]), and undoubtedly certain 
experimental efforts should be made to resolve this problem. 

Using our theoretical estimations for the decay properties of these nuclei (see above) we 
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have found that increase of the neutron fluence in the individual pulse by about three orders 
of magnitude as compared with existing pulsed reactors, i.e. up to 10^° netrons/cm^, could 
be quite sufficient to bypass the both gaps (see the solid curve in Fig. [7]). The details of 
the calculations and estimations of the yields of long-lived superheavy nuclei produced in 
neutron capture processes will be published elsewhere. 



III. SUMMARY 

We have shown for the first time that a macroscopic amount of the long-living superheavy 
nuclei located at the island of stability may be really produced in multiple (rather "soft") 
nuclear explosions, if such processes would be realized technically. This goal could be also 
reached by using the pulsed nuclear reactors of next generation, if their neutron fiuence per 
pulse will be increased by about three orders of magnitude. The experimental study of the 
decay properties of heavy nuclei located along the beta-stability line (and to the right of 
it, see Fig. [3]) is extremely important for a more accurate analysis of the neutron capture 



processes (including the astrophysical ones 20|]) in the upper part of the nuclear map. 



We are indebted to the DFG - RFBR collaboration and to GSI for support of our studies. 
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